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Abstract 

Many of the distinctive and subtle features of the dynamics in the £7a(1) channel in 
QCD can be related to gluon topology, more precisely to the topological susceptibility 
X{k 2 ) = ifd 4 x e ikx (0\T Q{x) Q(0)|0), where Q = ^trG^G^ is the gluon topological 
charge density. The link is the £^4(1) axial (ABJ) anomaly In this lecture, we describe 
the anomalous £^4(1) chiral Ward identities in a functional formalism and show how two 
apparently unrelated '£7^(1) problems' - the mass of the 7/ and the violation of the Ellis- 
Jaffe sum rule in polarised deep-inelastic scattering - can be explained in terms of the gluon 
topological susceptibility. They are related through a Ua(1) extension of the Goldberger- 
Treiman formula, which is derived here for QCD with both massless and massive quarks. 
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1. Introduction 

The central theme of this school is symmetry, in particular the realisation of 'hidden' 
or 'spontaneously broken' symmetries in quantum field theory. Another especially inter- 
esting realisation of symmetry in QFT is exemplified by the flavour singlet chiral Ua(X) 
symmetry in QCD with massless quarks. This is a subtle case, because although global 
chiral Ua(X) 1S a symmetry of the classical lagrangian of QCD, the corresponding current 
is not conserved in the quantum theory due to the well-known ABJ anomaly [1]. Never- 
theless, this 'anomalous symmetry' has important implications for the phenomenology of 
QCD. 

In this lecture, we describe how to analyse the consequences of the {7a (1) symmetry 
in QCD, using the formalism of anomalous chiral Ward identities. This formalism is then 
used to relate two '{7^(1) problems' - the mass of the rj' and the violation of the Ellis-Jaffe 
sum rule [2] in polarised deep-inelastic scattering (DIS) - to gluon topology in the form of 
the gluon topological susceptibility 

X (k 2 ) = tjd 4 x e ikx (0\T Q(x) Q(0)|0> (1.1) 

where Q = ^trG^G^ 1 ' is the gluon topological charge density. 

The topological susceptibility is a fundamental correlation function in pure gluody- 
namics (Yang-Mills theory) or QCD itself and is the key to understanding much of the 
distinctive dynamics in the Ua(X) channel. It has been studied non-perturbatively using 
lattice gauge theory [3-5], spectral sum rules [6-8], instanton models of the vacuum [9], etc., 
with the following results: 

Gluodynamics X FM (0) ~ -(180 MeV) 4 Xym(°) - -( 10 MeV) 2 

QCD (m q = 0) x(0) = x'(0) ~ (26 MeV) 2 (1.2) 

where x'(0) = ~^2X{k 2 )\ k _ - In gluodynamics, the value x YM (0) is well-established [3, 6] 
and can be calculated in lattice gauge theory. In QCD with massless quarks, the result 
x(0) = is an exact identity, following (as we show below) directly from the anomalous 
chiral Ward identities. The quoted value for x'(0) in full QCD was obtained using spectral 
sum rules[7,8]. Lattice techniques are still being refined[4,5] to produce a reliable result. 
The connection with topology arises as follows. Q is a total divergence, viz. 

Q = d^K^ (1.3) 

where = ^e^ pr7 tr(A L 'G pa - \gA u [A p , A a ]) is the Chern-Simons current. However, 
the integral over (Euclidean) spacetime of Q need not vanish. In fact, for gauge field 
configurations which become pure gauge at infinity, 

= „ 6 Z (1.4) 
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where the integer n is the topological winding number (technically, an element of the 
homotopy group n^(SU(N c )), where SU(N C ) is the gauge group[10, ch. 23]) or 'instanton 
number'. Instantons are classical field configurations which contribute to the path integral 
for which eq.(1.4) gives n^O. 

The connection with the UaO-) current in QCD arises through the famous ABJ 
anomaly. In a sense we shall make precise shortly, the flavour singlet axial current J^ 5 is 
not conserved in the quantum theory, even though it is the Noether current for the classical 
QCD lagrangian with massless quarks, but satisfies the anomaly equation (for m q = 0) 

d»J% - 2n f Q ~ (1.5) 

where n/ is the number of quark flavours. This identity provides the link between the 
quark dynamics and gluon topology. 

Our first example of a 'C/a(1) problem' concerns the mass of the rj'. This is much 
larger than would be the case if it were the pseudo-Goldstone boson for spontaneously 
broken Ua(X), as we would expect in the absence of the anomaly. Indeed, for massless 
QCD, the r]' would be an exact, massless Goldstone boson but for the anomaly. In section 
3, we explain how the Goldstone theorem is circumvented by the anomaly and derive the 
Witten-Veneziano mass formula [11, 12] for the r\' . This relates, to leading order in the 1/N C 
expansion, the mass of the r\' to the topological susceptibility in pure gluodynamics.W For 
QCD with massless quarks, the formula is simply 

^6 



77V - ^\/x YM (0) (1-6) 



The second example concerns the much-publicised 'proton spin' problem (for a review, 
see e.g. ref.[16]), i.e. the violation of the Ellis- Jaffe sum rule observed in measurements of 
the first moment of the polarised proton structure function g\ in deep inelastic scattering 
(DIS). As explained in section 4, the first moment J dxgf(x;Q 2 ) is related, through the 
OPE for the product of two electromagnetic currents, to the proton matrix elements of 
the SU(3) flavour singlet and non-singlet axial currents (p\J£ 5 \p), for a = 0,3,8. The 
corresponding 'axial charges' of the proton are denoted by a°(Q 2 ), a 3 and a 8 , where the 
singlet charge (only) has an explicit dependence on the scale Q 2 of the deep inelastic 
process. This is due to the non-trivial renormalisation of the singlet axial current induced 
by the anomaly. We have found[18,19] the following formula relating the ratio of the flavour 
singlet and non-singlet axial charges to the slope of the gluon topological susceptibility: 

The similarity to eq.(1.6) is striking. This formula is the basis of our quantitative resolution 
of the 'proton spin' problem[7], which is explained as a consequence of topological charge 
screening by the QCD vacuum. 



( 1 ) In this lecture, we derive this and other results directly from the anomalous chiral Ward identities. 
For the related approach of using effective chiral lagrangians in the 1/N C expansion to study rj' physics, 
see e.g. refs.[13,14,15]. 
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Although apparently quite different, these two examples are in fact intimately related 
through a U A (1) extension of the familiar Goldberger-Treiman relation. The essence of the 
conventional GT relation is that it links the dynamics of the pseudovector and pseudoscalar 
channels. This is just what we need to explain the OZI violation at the heart of the 'proton 
spin' problem, provided an extension of the GT relation to the anomalous Ua(1) can be 
found. This was achieved in refs. [17,18,19] and underlies our result (1.7). In section 5, we 
derive the Ua(1) GT relation. For massless QCD, we find 

2m N G° A = 2n fV ^(oj r v ONN\ k=0 (i-8) 

(Here, G\ is just an alternative notation for the nucleon axial charges, with a more con- 
venient SU(3) normalisation.) T v onn is a vertex describing the coupling of the nucleon 
to the rp - an unphysical state in QCD which in the OZI limit (see section 3) is identified 
with the exact UaO-) Goldstone boson. For QCD with massive quarks and flavour SU(3) 
breaking, the corresponding formula mixes the flavour singlet and non-singlet sectors. In 
this case, we find [8] 

2m N G A = F ab T v b NN \ k=Q (1.9) 

where F is determined from 

F ac Fl = \im-^ijdx e ik * (0|T d*Jfc(x) ^J* 5 (0)|0> (1.10) 

providing a natural generalisation of eq.(1.8). 



2. Chiral Ward Identities and the Renormalisation Group 

To study the phenomenology of the Ua(X) anomaly, we will need to find the Green 
functions of the composite operators which couple to the relevant physical states. These 
are the currents and pseudoscalar operators J^ 5 , Q, 0§ together with 4> a where 

(y ~ 

Jfc = W^q Q = ^trG M ,G^ 

<f> a 5 =q l5 T a q r = qT a q (2.1) 

is the field strength for the gluon field A M . In this notation, T % = |A* are flavour 
SU(nf) generators, and we include the singlet UaO-) generator T° = 1 and let the index 
a = We will only need to consider fields where i corresponds to a generator in 

the Cartan sub-algebra, so that a = 0,3, 8 for nj = 3 quark flavours. We define d- 
symbols by {T a ,T b } = d a b c T c . Since this includes the flavour singlet Ua(X) generator, 
they are only symmetric on the first two indices. For rif = 3, the explicit values are 
dooo = d 33 = d 88 = 2*^330 = ^880 = 1/3,(^338 = ^383 = -^888 = l/v 7 ^. (For further 
notation and description of the formalism used here, see refs. [19, 8].) 
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The Green functions, or correlation functions, are constructed from the generating 
functional VF[V^ 5 , 9, S§, S a ], where V® 5 , 9, S§, S a are the sources for the composite opera- 
tors J^5, Q, 05, (f) a respectively. Functional derivatives of W yield Green functions which 
are '1PI' w.r.t. the designated fields (composite operators). Explicitly, 



= J VAVqVq expz J d A x(C QC v + VfrJfc + 9Q + S 5 >g + S a cf> a ) (2.2) 



The chiral Ward identities are found by exploiting the invariance of W under a change 
of variables in the path integral corresponding to a chiral transformation q — > e taTai5 q. 
This gives 

J VAVqVq |>J»5 - 2n f 5 a0 Q - d adc m d (j) c 5 ~ s(i J g^qcd)] exp« J d A x(. . .) = 

(2.3) 

The terms in the square bracket are simply those arising from Noether's theorem, which 
defines the symmetry current through a functional derivative of the action, with the addi- 
tion of the anomaly term. This can be understood as arising from the non-invariance of the 
path integral measure VqDq in a background gauge field A^. A careful derivation of the 
anomaly by this method can be found in several standard textbooks on QFT, e.g. ref.[10, 
ch. 22]. The chiral variation term w.r.t. the elementary fields is then simply re-expressed as 
a variation w.r.t. the sources, giving finally the functional form of the (anomalous) chiral 
Ward identities: 

Q SW n . 5W , d 5W , ^dW , pd 5W n 

3 »W\ ~ 2nf6a0 l9~ ~ dadcm + dadcS W c ~ dadcS *Js~c = ( 2 - 4 ) 

This is the key equation which will be the basis of all the results derived in this lecture. 
It makes precise the anomaly equation (1.5). 

It will be useful in what follows to introduce some streamlined notation. The quark 
mass matrix is written as m a T a , so that for n/ = 3, 



m, 









m d 








m s 



(2.5) 



The chiral symmetry breaking condensates may be similarly written as 

(uu) \ 1 
(dd) = -<f>°l + 2(j) 3 T 3 + 2(j) 8 T 8 (2.6) 
(ss) J 3 

where ((p c ) is the VEV (qT c q). It is also convenient to introduce the still more compact 
notation 

M ab = d acb m c $ a6 = d abc {(j) c ) (2.7) 
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for which ^ 

-detM = m u md,m s -det<E» = (uu) (dd) (ss) (2-8) 

8 6 

Notice that the derivation of the chiral Ward identities sketched above was entirely in 
terms of the bare operators. Renormalised composite operators are defined as follows[20]: 

tO _ rv tO ja^O _ ja^O 

J ^5R — ZjJ fi5B j li5R ~ J n5B 

Q R = Q B -^-{l-Z)d»J% B 

<t>% R = Z^Ib <t>R = Z^ a B (2-9) 

where Z^ is the inverse of the mass renormalisation, = Z^ 1 . The anomalous dimensions 
associated with Z and Z^ are denoted 7 and 7^ respectively. Notice the mixing of the 
operator Q with d^J®^ under renormalisation. Most importantly, with these definitions 
the combination <9 M J^ 5 — 2n/Q occurring in the £7a(1) anomaly equation is RG invariant. 
The chiral Ward identities therefore take precisely the same form expressed in terms of the 
bare or renormalised operators. From now on, therefore, we use eq.(2.4) as an identity for 
renormalised composite operators (omitting the label 'R' for notational simplicity). 

It is also convenient to use a condensed notation where a functional derivative is 
represented simply by a subscript, with a spacetime integral assumed where appropriate. 
Also transforming to momentum space, we therefore write eq.(2.4) compactly as 

ikpWvfr - 2n f d a0 W 9 - M ac W Si + d adc S d W st - d adc SlW Sc = (2.10) 

The Ward identities for composite operator Green functions are derived by taking func- 
tional derivatives of this basic identity. We will need the following identities for 2-point 
functions: 

ikpWyayb ~ 2nfS a0 W e yb ~ M ac W S cyb^ = 

ikfjWy^e - 2nf5 a0 Wee - M ac W s ^e = 

ik^Wya S b - 2n f 5 a0 W 9S b - M ac WW - $ ab = (2.11) 

fio 5 o 5 5 

Combining the individual equations in (2.11), we find the important identity: 

k^Wya^ - M ac $ cb = Wga^ (2.12) 

where Sf) is the source for the current divergence operator D a = 2nf5 a0 Q + M ac (p^. In 
canonical notation, 

W Sf> S b D =ijdx e ik * (0|T D a (x) D b (0)|0> 

= ifdx e ikx (0|T d»j; b {x) ^J^ 5 (0)|0> (2.13) 



5 



The zero-momentum Ward identities play a special role. These follow immediately 
from eqs.(2.11) under the assumption that there are no massless particles (in particular, 
no exact Goldstone bosons) contributing 1/k 2 poles in the 2-point functions. With this 
assumption, we find simply 

2n f S a0 We9 + M ac W s ^e = 

2n f 5 a0 W dSb + M ac W S c S b + $ a6 = (2.14) 

5 5 5 

Combining these, we find that the topological susceptibility x(0) = Weo(0) satisfies the 
identity 

(2n f ) 2 X (0) = M ob M 0c W sbsc + M ob $ ob (2.15) 

5 5 

Another key ingredient in the discussion of the 'proton spin' problem and GT relations 
in sections 4 and 5 is the use of proper vertices for this set of operators. These are defined 
as functional derivatives of a generating functional T, which is itself constructed from W 
by a partial Legendre transform in which the transform is made only on the fields Q, (p^, <p a 
and not on the currents. The resulting proper vertices are 1PI w.r.t. the propagators for 
these composite operators only. As explained fully in refs. [19,21], by separating off the 
particle poles in the propagators, this is the definition which gives the closest identification 
of these field-theoretic vertices with physical low-energy couplings such as e.g. g^NN- We 
therefore define the generating functional r[V^ 5 , Q, 4> a ] as: 

r[v; 5 , q, c <p a ] = w\v^ e, s%, s a ] - J dx (oq + s^l + s a r) (2.16) 

The chiral Ward identities corresponding to eq.(2.10) are therefore: 

ik^Y V a s - 2n f 5 a0 Q - M ac (j) c 5 + d acd <j) d T^ - d ocd #?I> = (2.17) 

The Ward identities for the 2-point vertices will also be important. These follow directly 
from eq.(2.17): 

iku^ya yb + <§> ac F ±cyb = 

ikfjTy^Q - 2Uf5 a O + *o c r^Q = 

ik^Tya^ + ^acT m -M ab = (2.18) 

It is then straightforward to derive the following important identity, analogous to eq.(2.12): 

k^kjjYya yb + M ac $ cb = $ ac rV,d$ db (2.19) 

f-LO f5 5 o 

The renormalisation group equations for these quantities also play a key role in under- 
standing the physics of the {7a(1) channel. We therefore include here a brief and somewhat 
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novel discussion of the RGEs for the Green functions and proper vertices of these composite 
operators in a functional formalism. For further details, see especially refs. [22,19]. 

The fundamental RGE for the generating functional W follows immediately from the 
definitions (2.9) of the renormalised composite operators. It is: 

VW = -y{v% - ^- f 9^)w v o B + ^(s^Ws* + S a W s «) +... (2.20) 

where V = (v>-§ji + @-§jj — 7m Z~2 q m q~&k~) • ^ ne notation + . . . refers to the addi- 

tional terms which are required to produce the contact term contributions to the RGEs for 
n-point Green functions of composite operators. These are discussed fully in refs. [22,19], 
but will be omitted here for simplicity. They vanish at zero- momentum. 

The RGEs for Green functions are found simply by differentiating eq.(2.20) w.r.t. the 
sources. Simplifying the results using the chiral Ward identities (2.11), we find a complete 
set of RGEs for the 2-point functions. These are: 

VWyO yO = 2 r )WyO yO + . . . VWyO yb = 'jWyO yb + . . . VWya yb = + . . . 

VWyo e = 2-yWyo e + 7—*— M ob W v o s b + ■■■ 
VW v « 5 e = Wv^e + -y^-MobWyo^ +... 

VWyO^b = ( 7 + l4>)Wy^ Sl + . . . VWy^s* = l^Wy^b + . . . 

VW ee = 2*yW ee + 2 l7 —M 0b W eS b + ... 

2nf 

VW eS b = ( 7 + ^)W eS b + 7^- (M 0c W s = S 6 + $ 0& ) + • • • 

VW S a S b = 2^W s « S b + ... (2.21) 

It is straightforward to check the self-consistency of these RGEs with the Ward identities 
(2.11) and (2.14). The pattern of cancellations which ensures this is nevertheless quite 
intricate. 

Next, we need the RGE for the generating functional of the 1PI vertices. This follows 
immediately from its definition in eq.(2.16) and the RGE (2.20) for W: 

= 7 (v„° 5 - ^-r a„) rv„ s - lt (#r*. + <ri>) + . . . (2.22) 
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The RGEs for the 1PI vertices are found by differentiation, and using the Ward 
identities (2.18) to simplify the results, we find for the pseudoscalar sector: 



VT QQ = -2 7 r QQ + 2 7 



1 



2n f 



VT Q^ = -(7 + 7*)r ^ +7 

vv <t>«4>\ = - 2 7<^W +7 



QocTqqT^q 
1 



+ 



1 



2n f 



^Oc^QQT^b + TQ^TQ^b) - M ob T Q Q 
^oJVgQr^b - MobTQQ + a <-> 6 j + . . . 



+ . . . 

(2.23) 



Here, X> = X> + 7^>(</> a )^-. As explained in ref.[19], this is identical to the RG operator T> 
defined above (acting on W) when the sources are set to zero and the fields to their VEVs. 

It will also be useful to know the RGEs for the 3-point vertices coupling a pseudoscalar 
operator to the nucleon. In the same way, we find[19,8] 



OTqjvjv = ~1^qnn + 7 



2n f 



T ) T ( t > a NN — — 7</,rVjvjv + 7 



2n f 



®ob (r Q(f> bfQ NN + TQ Q t^b NN ^ + . .. 

^Ob^r^a^bf qjvJV + ^<t>tQ^<t> b 5 NNj - M 0a fQ NN 



+ ... 

(2.24) 



These RGEs play two roles in the discussion that follows. First, they will be used as 
consistency checks on the various formulae we derive. Second, and most important, they 
will provide the clue to identifying quantities which are likely to show violations of the OZI 
rule and those for which we may reasonably expect the OZI limit to be a good approxi- 
mation. This is because we can identify quantities which will be particularly sensitive to 
the Ua(X) anomaly as those which have RGEs involving the anomalous dimension 7. 



3. Pseudoscalar Mesons and the Witten-Veneziano Mass Formula for ?/ 

The anomalous chiral Ward identities have some immediate and important conse- 
quences. For simplicity, we specialise to massless QCD in this section. Then, integrating 
eq.(2.4) and evaluating with the sources set to their physical values (i.e. zero, apart from 
the source 9(x) which becomes the QCD theta angle #), we find 

SW ^#1 = (3.1) 



dd J 59{x) 

That is, massless QCD is independent of the theta angle. In fact [10, ch. 23], the same 
conclusion holds if any of the quark masses were to vanish. Then, the theta angle would 
have no effect and the strong CP problem would be automatically resolved in QCD. This 
is, however, an unrealistic solution since even m u 7^ 0. 
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It is equally simple to establish that the (zero-momentum) topological susceptibility 
vanishes for massless QCD. This can be read off immediately from the M = limit of 
eq.(2.15). There is, however, one subtlety here which is worth noticing. Writing the 
second identity in eq.(2.11) in canonical form, we have 

ifc"(0|T J° 5 Q\0) - 2n f (0\T Q Q\0) = (3.2) 

If there is no massless pseudoscalar meson (e.g. a Ua(X) Goldstone boson) coupling to the 
current, then clearly the first term vanishes at zero momentum. However, in this case the 
same conclusion follows even if there does exist such a particle, since the anomaly equation 
implies that the coupling of Q to this massless boson would vanish on-shell. In either case, 
we deduce 

X (0) = (0|TQQ|0)| fc=0 = (3.3) 

The next question is what happens to Goldstone's theorem in the presence of the 
anomaly. For the non-anomalous flavour non-singlet currents, the third identity in eq.(2.11) 
shows as usual that (for M = 0) if there is a symmetry breaking VEV $ ^ then there 
must exist, by Goldstone's theorem, a massless boson coupling derivatively to the current. 
What about the Ua(X) current? The relevant Ward identity, in canonical form, reads 

ik»(0\T J% 0j>|O> - 2n f (0\T Q <p° 5 \0) = 2(0°) (3.4) 

Since the r.h.s. is non-zero, in the absence of the anomaly term the only way the identity 
can be satisfied at k = is if there exists a massless Goldstone boson coupling to the 
current. However, the presence of the extra correlation function involving Q means that 
this conclusion no longer holds. The Goldstone theorem is evaded by virtue of the anomaly 
and there is no physical massless Ua(1) Goldstone boson. 

At its most basic, this is the resolution of the famous C/a(1) problem. However, things 
are of course not so simple. Recall from eq.(1.3) that Q may be written as the divergence 
of the (gauge non-invariant) Chern-Simons current K^. We can therefore construct a 
conserved, but gauge non-invariant, current J° 5 as follows 

J° 5 = J% - 2n f K fl (3.5) 

which satisfies the chiral Ward identity 

z^(0|T J° 5 0°|O) = 2(0°) (3.6) 

Applying Goldstone's theorem naively, we would then deduce that after all there must 
exist a massless boson in the QCD spectrum. However, this conclusion is false, although 
the precise reasons are still the subject of some debate. A full technical analysis of this 
Ua(X) problem would need (at least) another lecture, so instead we simply refer to the 
literature[23,24]. In brief, however, there are two possible escape routes (both of which 
may be true in a sufficiently precise formulation of the problem): 
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(i) the boundary conditions at spatial infinity on the behaviour of the gauge- variant current 
J° 5 imposed by the vacuum structure of QCD allow eq.(3.6) to be satisfied without coupling 
to a massless boson; 

(ii) a massless boson does indeed exist; however, it decouples from the positive-norm 
Hilbert space through a mechanism known as 'quartet decoupling'. This is essentially the 
same mechanism as is responsible for the decoupling from the physical spectrum of the 
ghosts and longitudinal and scalar components of the photon field in QED formulated in 
a covariant gauge. Its possible application to the Ua(X) problem has been elaborated by 
Kugo[25], who shows how the required Goldstone quartet is constructed by acting with 
the BRS operator Qb on the gauge- variant Goldstone field coupling to J° 5 . It represents 
a more realistic generalisation of the popular Kogut-Susskind dipole hypothesis [26] and 
reduces to it in the one case where the Kogut-Susskind hypothesis has been proved to 
work - the 2 dim Schwinger model. In this model, the 'quartet' splits into two 'dipoles' 
which independently cancel from the physical spectrum. 

Having convinced ourselves there is no real paradox associated with the Goldstone 
theorem applied to the anomalous {7a(1) current, we can ask whether it is possible to do 
better and identify the mass of the 'would-be Goldstone boson' rj' in some way with the 
anomaly. The answer is provided, in the context of the 1/N C approximation to QCD, by 
the Witten-Veneziano mass formula[ll,12]. 

We know that in nature the lightest pseudoscalar meson in the flavour singlet channel 
coupling to the current J^ 5 is the rj'. It follows from the second Ward identity in eq.(2.11) 
that this must also couple to the operator Q, producing a pole at k 2 = m 2 , in the topo- 
logical susceptibility Wqq. However, we have also seen that Wqq vanishes at k = 0. A 
reasonable parametrisation of the correlation function is therefore 

w >»^ = (^F^W^ 2 ' (3 ' 7) 

where A(k 2 ) is a pole-free, and therefore relatively smooth, function below the glueball 
threshold. A(k 2 ) is otherwise unconstrained by the chiral Ward identities. 
The residue of the pole evidently satisfies 

\(mv')\ 2 = j^ni^Aim 2 ,) (3.8) 

and so defining the RG non-invariant 'decay constant' f v ' by (see refs. [19,21] for a careful 
discussion of this point) 

(Q\4M = ik ^V ( 3 - 9 ) 
and using the anomaly equation, we find 

f 2 m%, = Aim*,) (3.10) 

Notice that neither side is a RG invariant, each scaling with the anomalous dimension 27. 
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We now introduce the large iV c approximation. To leading order in 1/N C , the anomaly 
is absent and the rj' mass vanishes. Formally, this is implemented by assuming that m v i is 
0(1/N C ). In fact, this can be alternatively re-formulated by using the OZI limit, i.e. the 
approximation to QCD in which the OZI rule becomes exact. The OZI limit is precisely 
defined [27] as the truncation of full QCD in which non-planar and quark-loop diagrams are 
retained, but diagrams in which the external currents are attached to distinct quark loops, 
so that there are purely gluonic intermediate states, are omitted. (This last fact makes 
the connection with the familiar phenomenological form of the OZI, or Zweig, rule.) This 
is a more accurate approximation to full QCD than either the leading large 1/N C limit, 
the quenched approximation (small n/ at fixed iV c ) or the leading topological expansion 
(iV c — > oo at fixed nf/N c . In the OZI limit, the C/a(1) anomaly is absent, as is meson- 
glueball mixing, and there is an extra {7a (1) Goldstone boson. 

Applying the large N c limit to the l.h.s. of eq.(3.7), we therefore find 

lim lim W ee (k 2 ) = j- 1 , o A(0) (3.11) 
fc ^0 7V c ^oo v ' (2n/) 2 v ' v ' 

(Notice that the zero- momentum and large N c limits do not commute.) We now put 
eqs.(3.10) and (3.11) together. To leading order in 1/N C , we may write f v > ~ y/2nff 7T 
and A(m v t) ~ A(Q). On the other hand, we may identify the l.h.s. of eq.(3.11) as the 
topological susceptibility in pure gluodynamics, since as noted above, mixing between the 
meson and glueball sectors vanishes to leading order in 1/N C , or equivalently in the OZI 
limit. 

Putting this together, we finally deduce the Witten-Veneziano formula for the mass 
of the 77', valid to leading order in the 1/N C expansion: 

J IT 

As discussed in the introduction, this is our first example of an explicit relation linking 
quark dynamics in the {7a (1) channel to gluon topology in the form of the topological 
susceptibility of pure gluodynamics. 



4. Topological Charge Screening and the Ellis-Jaffe Sum Rule 

The 'proton spin' problem, i.e. the question of why the first moment of the flavour 
singlet component of the polarised proton structure function g\ is anomalously suppressed, 
has inspired an impressive research effort, both theoretical and experimental, for over a 
decade. A recent review of this whole topic from the viewpoint adopted here can be found 
in ref.[16]. As is well-known from standard DIS theory, the first moment of g\ can be 
expressed in terms of the axial charges of the proton as follows: 

r? = f dx &{x;Q 2 ) = ^ s (a s ) (a 3 + ^a 8 ) + icf(a fl ) a°(Q 2 ) (4.1) 
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Here, Cf s , Cf are the appropriate Wilson coefficients arising from the OPE for two 
electromagnetic currents, while the axial charges are defined from the forward matrix 
elements of the axial currents with the normalisations 

{p, s\J* 5 \p, s) = (p,s\J 8 5 \p,s) = (p, s\ J%\p, s) = a°(Q 2 )s M (4.2) 

where s M = u(p, s)7 M 7514(7?, s) is the proton polarisation vector, a 3 and a 8 are known in 
terms of the F and D coefficients from beta and hyperon decay, so that an experimental 
determination of the first moment of g\ in polarised DIS allows a determination of the 
singlet axial charge a°(Q 2 ). The 'proton spin' problem is the fact that it is found exper- 
imentally that a°(Q 2 ) is strongly suppressed relative to a 8 , which would be its expected 
value if the OZI rule were exact in this channel. 

DIS is normally described theoretically using the QCD parton model. In this model, 
the axial charges are represented (in the AB renormalisation scheme [2 8, 2 9]) in terms of 
moments of parton distributions as follows [28]: 

a 3 = Au-Ad a 8 = Au + Ad-2As a°(Q 2 ) = Au + Ad+ As-n f ^Ag(Q 2 ) (4.3) 

In the parton model, the 'proton spin' problem takes the following form. In the naive, or 
valence quark, parton model we would expect the strange quark and gluon distributions to 
vanish, i.e. As = 0, Ag(Q 2 ) = 0. In that case, a = a 8 , the OZI prediction. Inserted into 
eq.(4.1), this gives the Ellis-Jaffe sum rule [2]. However, the observed suppression a°(Q 2 ) < 
a 8 can be accommodated in the full QCD parton model by invoking either or both a non- 
zero polarised strange quark distribution As 7^ or a non-zero polarised gluon distribution 
Ag(Q 2 ) 7^ 0. An interesting conjecture (in line with the insights of the approach discussed 
here) is that the suppression is primarily due to the gluon distribution [28], although a 
quantitative prediction would still only follow if Ag(Q 2 ) can be independently measured, 
either through a precise analysis of the Q 2 dependence of g\ [29] or directly through other 
less inclusive high energy processes such as open charm production. Notice, however, 
that even in the QCD parton model picture, it is not possible to identify a°(Q 2 ) with 
spin[17,16]. This identification only holds for free quarks, in which case the Q 2 scale 
dependence, which is related through eq.(2.9) (see also eq.(5.26)) to the Ua{X) anomaly, 
disappears from a . It must be emphasised that the so-called 'proton spin' problem is not 
a problem of spin - rather, it is a question of understanding the dynamical origin of the 
OZI violation a°(Q 2 ) <a 8 . 

In this section, we shall discuss a less conventional approach (the 'CPV method[16]) 
to DIS based on the composite operator propagator - proper vertex formalism described 
in section 2. The starting point, as indicated above, is the use of the OPE in the proton 
matrix element of two currents. This gives the standard form for a generic structure 
function moment: 

f dxx n - 1 F{x;Q 2 ) = Y,CnQ 2 )(p\O r tmp) (4.4) 
Jo 
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where Of are the set of lowest twist, spin n operators in the OPE and C™(Q 2 ) the corre- 
sponding Wilson coefficients. In the CPV approach, we now factorise the matrix element 
into the product of composite operator propagators and vertex functions, as illustrated in 
Fig. 1. 




Fig.l CPV description of DIS. The double line denotes the composite operator propogator and the lower 
blob the 1PI vertex. 

To do this, we first select a set of composite operators Oi appropriate to the physical 
situation and define vertices ^o iPP as 1PI with respect to this set. Technically, this is 
achieved as in eq.(2.16) by introducing sources for these operators in the QCD generating 
functional, then performing a (partial) Legendre transform to obtain a generating func- 
tional r[Cy . The 1PI vertices are the functional derivatives of . The generic structure 
function sum rule (4.4) then takes the form 




dxx™- 1 f(x,q 2 ) = ^^c] n) (g 2 )(o|T of a,|o)r a 



=EEw (4-5) 

i 3 

in a symbolic notation. 

This decomposition splits the structure function into three pieces - first, the Wilson 
coefficients cf(Q 2 ) which control the Q 2 dependence and can be calculated in pertur- 
bative QCD; second, non-perturbative but target-independent QCD correlation functions 
(0\TOf Oi\0); and third, non-perturbative, target-dependent vertex functions ^Q. pp de- 
scribing the coupling of the target proton to the composite operators of interest. The vertex 
functions cannot be calculated directly from first principles. They encode the information 
on the nature of the proton state and play an analogous role to the parton distributions 
in the more conventional parton picture. 

It is important to recognise that this decomposition of the matrix elements into prod- 
ucts of propagators and proper vertices is exact, independent of the choice of the set of 
operators Oi. In particular, it is not necessary for 0\ to be in any sense a complete set. 
All that happens if a different choice is made is that the vertices , pp themselves change, 
becoming 1PI with respect to a different set of composite fields. Of course, while any set of 
Oi may be chosen, some will be more convenient than others. Clearly, the set of operators 
should be as small as possible while still capturing the essential physics (i.e. they should 
encompass the relevant degrees of freedom) and indeed a good choice can result in vertices 
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^OiPP which are both RG invariant and closely related to low energy physical couplings, 
such as g-KNN- In this case, eq.(4.5) provides a rigorous relation between high Q 2 DIS and 
low-energy meson-nucleon scattering. 

For the first moment sum rule for g\ [19,7,16], it is most convenient to use the U a(S) 
anomaly equation immediately to re-express a°(Q 2 ) in terms of the forward matrix element 
of the topological charge Q, i.e. 



a°(Q 2 ) = 



2mN 



2n f {p\Q\p) 



(4.6) 



where mjv is the nucleon mass. 

Our set of operators Oi is then chosen to be the renormalised flavour singlet pseu- 
doscalars Q and $5, where $5 is simply the operator 0° of eqs.(2.1) and (2.9) with a special, 
and crucial, normalisation. The normalisation factor is chosen such that in the absence of 
the anomaly, or more precisely in the OZI limit of QCD (see section 3), $5 would have the 
correct normalisation to couple with unit decay constant to the C/a(1) Goldstone boson 
which would exist in this limit. This also ensures that the vertex is RG scale independent, 
as we prove in section 5 where we relate this discussion to the UaO-) Goldberger-Treiman 
relation. The vertices are defined from the generating functional (2.16). We then have 



1 singlet 



1 1 



9 2tm tv 



2n f Cf{a s ) 



(0|TQQ|0)f Qpp +(0|TQ$ 5 |0)f $ 



SPP 



(4.7) 



where the propagators are at zero momentum and the vertices are 1PI wrt Q and $5 only. 
For simplicity, we have also introduced the notation iuYQ pp u = Yq pp u^u, etc. This is 
illustrated in Fig. 2. 



XQ 



Q 



X Q 




<*>5 




Fig. 2 CPV decomposition of the matrix element (p\Q\p). 

The composite operator propagator in the first term is simply the (zero-momentum) 
QCD topological susceptibility x(0) which, as we have seen in sections 2 and 3, vanishes 
for QCD with massless quarks. Furthermore, with the normalisation specified above for 
$5, the propagator (0|T Q $s|0) at zero momentum is simply the square root of the slope 
of the topological susceptibility. 

To see this, notice that by virtue of their definition in terms of the generating functional 
(2.16), the matrix of 2-point vertices in the pseudoscalar sector is simply the inverse of the 
corresponding matrix of pseudoscalar propagators, i.e. 



QQ 
02Q 



r, 



Wee W ( 



W, 



9S° 

s° 5 s° 



(4.8) 
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This implies 

r m = -W ee {det Wss)' 1 (4.9) 

letting S represent the set {6>, S®}. Differentiating w.r.t. k 2 and taking the limit k 2 = 0, 
and exploiting the fact that Wee(0) vanishes, we find 

<l T ^\ k =o = x(0)W d - s l(0) (4.10) 



dk 2 



Finally, normalising the field $5 proportional to (fi® such that ^2-r^> 5 ^> 5 1 fc=Q = 1 (see also 
the discussion in section 5, following eq.(5.7)), we find the required relation 



\TQ^\0)\ k=Q = V7W) (4.11) 



We therefore find: 



1 1 



r? singlet = Q^ 2n f C lM VTW) f* 5PP (4.12) 

The slope of the topological susceptibility x'(0) is not RG invariant but, as shown in 
eq.(2.21), scales with the anomalous dimension 27, i.e. 

f t V7W) = 7 (4.i3) 

On the other hand, the proper vertex has been chosen specifically so as to be RG invariant, 
as proved in section 5. The renormalisation group properties of this decomposition are 
crucial to our proposed resolution of the 'proton spin' problem. 

Our proposal[19] is that we should expect the source of OZI violations to lie in the 
RG non-invariant, and therefore anomaly-sensitive, terms, i.e. in x'(0) rather than in the 
RG invariant vertex. Notice that we are using RG non-invariance, i.e. dependence on the 
anomalous dimension 7, merely as an indicator of which quantities are sensitive to the 
anomaly and therefore likely to show OZI violations. Since the anomalous suppression in 
is thus assigned to the composite operator propagator rather than the proper vertex, 
the suppression is a target independent property of QCD related to the anomaly, not a 
special property of the proton structure. 

To convert this into a quantitative prediction we use the OZI approximation for the 
vertex r$ 5Pp . In terms of a similarly normalised octet field $|, this is r$ 5PP 
Notice that the normalisation is crucial in allowing the use of the OZI relation here. 
The corresponding OZI prediction for a/V(0) would be fn/VQ- These OZI values are 
determined by comparing the result (4.12) (at least that part relating to the proton matrix 
element) with the conventional Goldberger-Treiman relation for the flavour octet axial 
charge in the chiral limit (see section 5). This gives our formula 

a ^=£vm* (4.14) 
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for the flavour singlet axial charge. Incorporating this into the formula for the first moment 
of the polarised structure function, we find 

rising = lc?(a 8 )a*^Vm (4-15) 

Finally, substituting this into the expression for the complete first moment and using 
our spectral sum rule derivation of x'(0), which gives a suppression factor of around 0.6 
(see refs. [7,8,16] for more details), we obtain our prediction in the chiral limit 

a°(Q 2 = lOGeV 2 ) = 0.33 ± 0.05 

T P (Q 2 = lOGeV 2 ) = 0.144 ± 0.009 (4.16) 

This is to be compared with the OZI (Ellis-Jaffe) prediction a = 0.58 ± 0.03 and the 
current experimental data from the SMC collaboration[30]: 

r?(Q 2 = 10GeV 2 )| = / dxg p 1 (x;Q 2 ) = 0.141 ±0.012 (4.17) 

^ ' ' ^0.003 

The result for the entire first moment depends on how the extrapolation to the unmeasured 
small x region x < 0.003 is performed. This is still a controversial issue. Using a simple 
Regge fit, SMC find = 0.142 ± 0.017 from which they deduce a = 0.34 ± 0.17, while 
using a small x fit using perturbative QCD evolution of the parton distributions[31] they 
find rf = 0.130 ± 0.017 and a = 0.22 ± 0.17 (all at Q 2 = lOGeV 2 ). 

More recently, SMC have published an alternative analysis [32] of their data, this time 
quoting a slightly lower number for the integral over the measured region of x: 

r?(Q 2 = 10GeV 2 )L >Q 003) = f dxg p (x;Q 2 ) = 0.133 ±0.009 (4.18) 

' ^0.003 

Clearly it is premature to draw too strong a conclusion given the large errors on the 
experimental determinations of T p and a and the uncertainty over the small x extrapola- 
tion. Nevertheless, it is extremely encouraging that our prediction is firmly in the region 
favoured by the data and gives us confidence that our explanation of the 'proton spin' 
problem in terms of the topological susceptibility is correct. (For a more extensive discus- 
sion of the phenomenological aspects of this approach to the 'proton spin' problem, and a 
proposal to test our 'target-independence' conjecture in semi-inclusive DIS experiments at 
polarised HERA, see refs. [16,33,34]). 

Finally, we should draw attention to the basic dynamical mechanism responsible for 
the suppression of the 'proton spin'. What we have shown is that when a matrix element 
of the topological charge is measured, the QCD vacuum screens the topological charge 
through the zero or anomalously small values of the susceptibility x(0) and its slope x'(0) 
respectively (see Fig. 2). The mechanism is analogous to the screening of electric charge 
in QED. There, because of the gauge Ward identity, the screening is given entirely by the 
('target independent') dressing of the photon propagator by vacuum polarisation diagrams, 
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leading to the relation = eB^/Zs (with Z3 < 1) between the renormalised and bare 
charges, in direct analogy to eq.(4.14) above with the topological susceptibility playing the 
role of the photon propagator. 



5. The UaO) Goldberger-Treiman Formula 2 



In this final section, we present a unified derivation of the Goldberger-Treiman rela- 
tions (1.9) for the flavour singlet and non-singlet axial charges for QCD with non- vanishing 
quark masses in the functional composite operator propagator - vertex formalism [8]. The 
zero-mass limit of the formula for the singlet £^4(1) charge reduces to the expression al- 
ready used in section 4 (see eq.(4.12)) in the context of the 'proton spin' problem. Notice 
that these new GT relations are more general than the conventional formulae. Indeed, 
they are exact field-theoretic results in QCD. The familiar PCAC forms (in the flavour 
non-singlet channels) are obtained by approximating the 1PI vertices by the correspond- 
ing low-energy meson-nucleon couplings and by approximating the slopes of the current 
correlation functions (1.10) by decay constants. Away from the chiral limit, both these 
approximations assume pole dominance of the matrix elements and correlation functions 
by the pseudo-Goldstone bosons. 

The axial charges G A are defined as the form factors in the forward nucleon matrix 
elements of the axial currents, viz. 



G a A s 



(5.1) 



{p,s\J^\p, s ) 

(Compare with the normalisations in eq.(4.2).) 

To express this matrix element in terms of composite operator propagators and the 
associated 1PI vertices, we first introduce an interpolating field N and source S n for the 
nucleon. (Notice that this is purely a formal device - there is no dynamics implicit in this 
step.) The matrix element is then just the 3-point function Wv a s N s N with the external 
propagators amputated. This can be re-expressed in terms of the vertex functional V as 
follows (see refs.[19,8] for a derivation of the relevant formulae involving partial Legendre 
transforms) : 



(p, s \ j »5\p, s ) = U (P' S ) 



W 



w v°s N s N W, 



/lx5 



SnS? 



u(p, s) 



= u(p,s) 



Fv£ 5 NN + Wv" 5 G ^QNN + Wya R gb T^b 



p.. 5 



/j.5 5 



NN 



u{p,s) (5.2) 



Since the propagators on the r.h.s. vanish at zero momentum (this requires the absence of 
any 1/k 2 poles, which as we have seen in section 3 is assured by the Ua(X) anomaly and 
quark masses), we find simply 



2m^G a A wy^u = u 



MV; 5 Ar;v| fc=0 



u 



(5.3) 



The material in this section was not presented in the lecture in Zuoz. However, it is closely related 
to the topics described in sections 2-4 and should help to clarify the discussion given there. 
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The GT relations then follow immediately from the Ward identity (2.17) for T. Differen- 
tiating w.r.t. the nucleon fields, we find 

2m N G a A = $ ab f^ NN \ k=Q (5.4) 

where as before we define %u T^a NN u = T^nn ^75^, etc. 

A non-forward version of the GT relation, which is closer to the analysis in section 
4, can also be found from eq.(5.2) by using the Ward identities (2.11) for the propagators 
together with (2.17) for Tv^nn- This allows us to write, for all k, 

2m N G\{k 2 ) + k 2 G%{k 2 )=-{2n f 5 a0 Wee+M ac W s ^e)t QNN 

- (2n f S a0 W eS b + M ac W s . s ^jt^ NN (5.5) 

Gp(k 2 ) is the pseudoscalar form factor in the non-forward matrix element (p, s\J® 5 \p, s), 
and again has no 1/k 2 pole for the reasons given above. This expression clearly reduces to 
eq.(5.4) on using the zero-momentum Ward identities (2.14) for the propagators. In the 
chiral limit (M = 0), this can be compared with eq.(4.7). 

The remaining step to convert eq.(5.4) into the useful form of the GT relations is to 
normalise the field 4>% appropriately. Clearly, eq.(5.4) is independent of the normalisation. 
However, with a suitable choice, the vertices can be made both RG invariant and essentially 
identical to the physical Goldstone boson couplings g^NN etc. To achieve this, we define 
normalised fields 

rf = B ab cf> b 5 (5.6) 
where B is a constant matrix such that 3 

This condition ensures that the fields rf have unit coupling to the Goldstone bosons. 

The case of the singlet rf is of course special, since it is only after mixing with the 
topological field Q (and then flavour mixing) that it becomes the physical rj'. In fact, this 
is why it is most convenient to impose the normalisation condition as above on the matrix 
of 2-point vertices r„ a „b, which is the inverse of the pseudoscalar propagator matrix, since 
this most simply characterises the rf before mixing with Q. The corresponding state is 
the unphysical 'OZI Goldstone boson' introduced in ref.[17] and extensively discussed in 
refs. [19,21]. 

It can now be proved that the vertices T v a NN defined with the fields normalised 
according to eq.(5.7) are RG invariant. The proof is based on the functional form of the 
RGEs introduced in section 2 and is given at the end of this section. 

Re-expressing eq.(5.4) in terms of the properly normalised vertices, we have 

2m N G a A = §a C Bj b r v b NN \ k=0 (5.8) 

3 Applied to the flavour singlet field <j>® in the chiral limit, this is the same condition used to normalise 
the field $5 in section 4. 
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where B is to be determined from the 2-point vertex condition 



• V 4>Z4>t I fc=n - B ac -TT2 V n c r 1 d I k =n B db - B Ic B cb (5.9) 



^2^5*5^=0 ac V^V a I fe=0 

The straightforward approach to finding T^a^b is to write it as one component of the 
inverse of the propagator matrix Wss (with S = {6,S%}). This was the approach used 
in section 4 to relate T^o^o to the topological susceptibility Wgg. However, inverting this 

^5^5 

matrix in the m q ^ 0, multi-flavour case is cumbersome, so here we use an indirect but 
more elegant method. First, we use the Ward identity (2.19) to express T^b in terms of 

Y\ra T/b 



V\V\ 
(AO i/o 



k^KYya yb = $ ac T^d $ db - M ac $ cb (5.10) 



Then, using a general identity for partial Legendre transforms[19,8], we relate Y V a V b to 

/a 5 i/5 

the current propagator W ya V b : 

jU.5 i/5 

Fya yb = Wya yb — Wya S W7^- W T yb (5-H) 

where 5 and T represent {6>, S5 1 }. Combining with eqs.(2.12) and (2.19) finally yields 

$ ac F^d <S> db = W S a S b + kf,W v? s W$t W TV bk v (5.12) 

' 5 ' 5 D D H' ' i/5 

and so, in matrix notation, 

® BTm = W? Ws ° s °\k=o + ^faWv^s Ws± W TVv5 K)\ k=0 (5.13) 

The argument is almost complete. $£? T is precisely the combination we need for the GT 
relations (5.8), and is related by eq.(5.13) directly to the first moment of the correlation 
function (2.13) for the divergences of the currents. 

It remains to show that the final term in eq.(5.13) vanishes. The first and last fac- 
tors are of 0(/c 2 ), so this will contribute zero unless there is a 1/k 2 pole in the inverse 
pseudoscalar propagator matrix W^. As already mentioned, there are no 1/k 2 poles in 
the propagators themselves, so all we need show is that the determinant A(k 2 ) of this 
propagator matrix is non- vanishing at k = 0. This follows from the formula 

A(0) = Wee(0) (detM)" 1 det$ (5.14) 

since detM and det$ are non-zero (eqs.(2.8)) and Wqq is non- vanishing away from the 
chiral limit. An elegant proof of eq.(5.14), based on the zero-momentum Ward identities 
for Wst is as follows. Define 
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with A = detW. Then, 



detM A = detM detW 



W, 



9S\ 



M ar W, 



(5.16) 



Using the zero-momentum Ward identities (2.14) and, in the case of the a = row taking 
a linear combination with the first (6) row, the determinant simplifies, leaving 



detM A 



Wee W t 

-$ab 



es b 5 



= Wee det$ 



(5.17) 



as required. 

This completes the proof of the GT relations. To summarise, we have shown that the 
flavour singlet and non-singlet axial charges are given by the single, unified relation 

2m^G\ = F ab T v b NN \ k=Q (5.18) 
where F = &B T is determined from 

ffT = Jr 2 W s»s D \ k=0 (5.19) 

that is, 

F ac Fl = Um i j dx e ikx (0|T d»J« 5 {x) ^4 5 (0)|0) (5.20) 

The current algebra approximation to eq.(5.20) is easily found. Consider just the G\ 
relation for simplicity, since in this case isospin invariance (assume m u = ma = 0) ensures 
that there is no flavour mixing. Applying pole dominance to the current correlator, we 
find 

* = £8s*(**f) = * (b ' 21) 

where in the intermediate step we have distinguished the factors of k 2 arising from the 
decay constant and pole terms. Under the same PCAC assumptions, the vertex can be 
identified with the low-energy pion-nucleon coupling constant, i.e. T v 3 NN\k=o = 9-kNN- 
Writing G\ in standard notation as where qa is the isotriplet axial- vector coupling 

of the nucleon measured in beta decay, we therefore recover in this approximation the 
standard GT formula 

fi^gnNN = m N g A (5.22) 

We emphasise again, however, that the relations (5.18) and (5.20) are exact - they are nei- 
ther dependent on this interpretation nor make any PCAC assumption or approximation. 
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The GT relation simplifies in the chiral limit, where flavour mixing is absent. In this 
case, the singlet axial charge is simply given by 

2m N G° A = 2n f y^O) f v o NN \ k=Q (5.23) 

where as usual x'(0) is the slope of the topological susceptibility. 

This formula can now be recognised as the basis of our resolution of the 'proton spin' 
problem (see eq.(4.12)). The general formula (5.18) extends this to include mixing with 
the flavour non-singlet sector, introducing a matrix structure which replaces the simple 
square root of x'(0) in eq.(5.23) and generalising the fields in the correlation function to 
the entire divergence of the current including mass terms as well as the anomaly Q. 

The final task is to determine the RG properties of the various quantities entering into 
the derivation of the GT relations, and verify the assertion that the 1PI vertices defined 
here are RG invariant. 

The RGE for the matrix B ab , which relates the <p^ fields to the canonically normalised 
yf fields by i] a = B ab (j)^ is readily found using the RGE (2.23) for the 2-point vertex 
Y^a^b. From the definition (5.9), using eq.(2.23) and the zero-momentum limit of the 
Ward identities (2.18), we deduce 

VB ab = - 1(j) B ab + 7 S ac $ c0 $o6 1 (5-24) 

The RGE for F ab now follows immediately from its definition F = <&B T and the RGE 
= 7<^$. It is simply 

VF ab = -f5 a0 F ob (5.25) 

The final step in proving RG consistency of the unified GT formulae is to show that 
the vertices T^nn (at k = 0) are RG invariant. Eq.(5.25) then ensures the required RGE 
for the axial charges, 

VG\ = -y6 a0 G a A (5.26) 

showing that only the singlet axial charge has a non-trivial RG scaling. To check this 
explicitly, notice that eq.(2.24) for T^nn simplifies at k = 0. The contact terms vanish 
and using the zero-momentum Ward identities (see eq.(2.18)) we find 

VY 4> a NN | k=Q = -7^r^ ArAr I fc=Q + 7$^ 1 $ 06 r^ A r A r I fc=Q (5.27) 

Since r<^g atat | ^._q = Bab^r) bNN \k=0' anc ^ com P arm g eq.(5.27) with the RGE (5.24) for £?, 
we confirm 

Vf v a NN \ k=Q = for all a (5.28) 
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